ChinaXiv 合 作 期 刊 


第 38 卷 第 6 期 干旱 区 研究 Vol.38 No.6 
2021 年 11 月 ARID ZONE RESEARCH Nov. 2021 


数值 模拟 方法 在 河流 冰 塞 水 位 模拟 中 的 应 用 
一 一 以 黄河 三 湖 河 口 弯 道 段 为 例 
RILA, BHR”, HAX, RREY, WAR, 谢 冰 绮 3 


(1. 河 海 大 学 水 文 水 资源 与 水 利 工 程 科 学 国家 重点 实验 室 ,江苏 南京 ”210098; 2. 河 海 大 学 水 文 水 资源 学 院 ， 
江苏 南京 ”210098; 3. 珠江 水 利 委员 会 珠江 水 利 科学 研究 院 , 广 东 广州 ”510630) 


摘 要 : 运用 开源 的 河 冰 数 值 模型 RIVICE 模拟 黄河 三 湖 河口 弯 道 段 2008 一 2010 年 开 河 期 的 冰 塞 水 位 ,并 进行 模 
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型 参数 的 敏感 性 分 析 。 研 究 所 使 用 的 流量 、 水 位 、 冰 厚 等 数据 均 来 自 水 文 年 鉴 。 结 果 表明 ,用 RIVICE 模 型 模拟 开 


a 


i 


期 冰 赛 水 位 具有 很 高 的 可 行 性 。 冰 塞 回 水 水 位 对 上 游 边界 流量 十 分 敏感 ,说 明了 水 利 工程 调控 的 重要 性 ,研究 


黄河 宁 蒙 河 段 由 于 特殊 的 地 理 位 置 水 力 条 
件 .河道 形态 以 及 人 为 因素 的 影响 ,是 凌 汛 灾害 频 
发 的 典型 河 段 HH。 凌 汛 的 发 生 是 由 于 冰 寨 ,国际 上 
将 冰 塞 划分 为 窗 河 冰 塞 和 宽 河 冰 塞 ,分别 是 水 力 
增 厚 和 机 械 增 厚 所 致 引 , 国 内 则 将 其 称 为 冰 塞 和 冰 
坝 ” 。 冬 季 冻 结 期 易 发 生 冰 塞 ,春季 开 河 期 易 发 生 
冰 坝 , 开 河 期 的 冰 坝 会 使 上 游 水 位 大 幅 雍 高 ,造成 
更 严重 的 淹没 和 演 堤 风险 。 据 统计 1991 一 2010 年 
黄河 宁 蒙 河 段 产生 的 局 部 冰 坝 凌 灾 多 达 36 场 , 单 次 
损失 最 多 达 7 亿 元 ”。2010 年 后 ,受到 上 游 水 利 工 
程 调控 的 影响 , 冰 坝 凌 灾 减少 ,多 以 文 开 河 为 主 呈 。 
三 湖 河口 弯 道 段 因 其 独特 的 河道 形态 ,属于 极 易 卡 
冰 位 置 ,虽然 在 上 游 水 利 工程 的 调控 和 人 工 破 冰 的 
作用 下 , 凌 汛 灾害 得 到 一 定 缓解 ,但 对 于 开 河 期 冰 
坝 洪 水 的 研究 需要 进一步 的 探究 。 

国外 对 冰 塞 的 机 理 研 究 始 于 20 世 纪 60 年 代 ， 
最 初 基于 稳 态 进行 冰 的 理论 研究 ”。 随 着 计算 机 技 
术 的 发 展 ,国际 上 有 很 多 学 者 建立 了 模拟 河 冰 的 数 
值 模型 ,如 一 维 模型 ICEJAMIG , RIVJAM" ,RiverID 5 , 
RIVICE ?' , RICEN 0 等 ,二 维 模型 DynalCE"" 和 
CRISSP2D'2。 国 内 在 对 河 冰原 型 观测 和 实验 人 研究 
等 方面 取得 了 很 多 成 果 天 5 ,但 对 河 冰 的 数值 模拟 
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通过 数值 模型 对 冰 塞 水 位 进行 模拟 及 演变 预测 ,为 水 利 工程 调控 和 闭 汛 防治 提供 一 定 的 参考 。 
关键 词 : 黄河 三 湖 河口 弯 道 段 ; 冰 塞 水 位 ; RIVICE 模型 


起 步 较 晚 ,研究 较 少 。 茅 泽 育 等 "结合 热力 学 和 动 
力学 计算 ,模拟 河道 结 冰 封 冻 过 程 中 冰 塞 体 的 演变 
以 及 对 水 位 的 影响 ; 素 鸿 兰 等 ”通过 耦合 水 温和 订 
盖 生 消 模 型 对 静水 状态 下 的 水 体 结 冰 融 化 过 程 进 
行 了 数值 模拟 ; EPUKA E RAZ V RICEN 模型 为 
基础 并 做 出 改进 ,将 使 冰 盖 破碎 的 临界 流量 由 1 个 
常数 改 为 1 个 经 验 公 式 , 开 发 了 综合 一 维 非 稳 态 河 
冰 动 力 模型 YRIDM; 李 超 " ”用 二 维 模型 DynaICE 模 
拟 了 三 湖 河口 弯 道 段 的 封 河 过 程 。 
不 同 于 已 有 研究 ,本 研究 的 目标 是 模拟 开 河 
,完整 冰 盖 破裂 之 后 向 下 游 流 去 ,在 河道 的 易 卡 
冰 位 置 堵塞 ,降低 河道 过 流 能 力 , 引 起 上 游 水 位 刻 
高 的 过 程 。 这 种 堵塞 可 能 是 由 单 层 浮 冰 构成 的 表 
HIKE” ,流量 大 时 ,也 有 可 能 受 力 增 厚 形成 冰 坝 
( 宽 河 冰 塞 ) ,所 以 ,本 文 将 这 一 过 程 统称 为 冰 塞 。 
本 研究 选取 RIVICE 模型 模拟 这 一 过 程 ,一 是 因为 
RIVICE 模 型 不 仅 可 以 模拟 冬季 河道 结 冰 封冻 ,也 可 
模拟 开 河 期 冰 盖 破碎 后 在 力 的 作用 下 增 厚 ,形成 阻 
水 冰 塞 "1; 二 是 因为 RIVICE 模 型 采用 冰 动 力 与 水 动 
力 双 向 松散 耦合 的 方式 来 模拟 冰 塞 ,不 需要 求解 复 
杂 的 微分 方程 ,计算 更 加 简单 高 效 2。 
本 研究 选取 黄河 内 蒙古 三 湖 河口 弯 道 段 为 研 
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究 河 段 ,模拟 开 河 时 期 因 河 道 卡 冰 而 造成 的 水 位 刻 
高 。 旨 在 验证 RIVICE 模型 模拟 的 可 行 性 ,以 及 为 
上 游 水 利 工程 在 开 河 期 的 防凌 汛 调 度 提 供 科 学 依 
据 ,在 尽 可 能 不 造成 严重 冰 塞 洪水 的 情况 下 ,向 下 
游 释 放 更 多 的 流量 。 


1 材料 与 方法 


1.1 研究 区 概况 

选取 地 处 半 干 旱 区 的 黄河 三 湖 河 口 弯 道 段 (图 
1) ,范围 是 亿 利 黄河 大 桥 以 上 的 9.89 km。 此 处 河床 
坡 降 小 ,河道 中 有 心 滩地 形 复杂 ,加 之 泥 沙 淤积 ,过 
流 能 力 弱 。 此 处 有 2 个 弯 道 、 弯 道 曲 率 大 ,第 2 个 弯 
道 处 河道 收缩 近乎 一 半 ,加 之 受 下 游 亿 利 黄河 大 桥 
桥墩 的 阻 水 影响 ,是 极 易 卡 冰 的 河 段 。 黄 河 三 湖 河 
口上 游 段 河流 呈 自 南 向 北 流向 ,因此 在 春季 开 河 
期 ,上 游 先 开 河 ,破裂 的 冰 盖 向 下 流 至 此 处 ,加剧 了 
形成 冰 塞 洪水 的 风险 ,所 以 在 此 段 进行 冰 塞 洪水 水 
位 的 预报 预警 研究 意义 突出 。 参 考 已 有 人 研究 ,该 
河 段 河 底 坡度 为 0.0001, 据 此 概 化 出 了 宽 浅 矩形 河 
道 断面 。 
1.2 数据 来 源 

本 研究 采取 2008 年 .2009 年 .2010 年 畅 流 期 和 
开 河 期 的 洪水 数据 来 率 定 和 验证 模型 ,洪水 数据 均 
来 自 水 文 年 鉴 。 模 拟 的 上 边界 设置 为 流量 过 程 , 下 
边界 设置 为 水 位 过 程 。 因 选取 河道 短 ,边界 处 没有 
水 位 和 流量 的 实测 数据 ,所 以 上 边界 的 流量 和 下 边 
界 的 水 位 都 用 三 湖 河口 站 的 数据 。 与 此 同时 ,用 三 
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湖 河口 站 的 数据 来 验证 模拟 结 
1.3 冰 塞 计算 理论 方法 

RIVICE 模型 是 由 加 拿 大 环境 与 气候 变化 委员 
会 (Environment and Climate Change Canada, ECCC) 
开发 的 开源 模型 。 该 模型 是 一 维 的 水 动力 与 冰 动 
力 过 程 的 双向 松散 耦合 模型 ,用 隐 式 有 限 差 分 格式 
对 圣 维 南 方程 组 进行 求解 ,并 加 入 了 河 冰 动力 过 
程 。 模 型 可 模拟 冰花 的 形成 . 岸 冰 的 演进 、 冰 盖 的 
形成 .来 冰 在 冰 羡 前 缘 的 并 置 销 冰 、 冰 的 输 移 、 冰 
盖 的 融化 与 破碎 、 冰 塞 和 悬 坝 这 些 河 冰 过 程 2。 

Lindenschmidt 等 ”已 将 此 模型 用 于 加 拿 大 的 
阿 萨 巴 斯 卡 河 .和 平河 呈 等 河流 的 冰 塞 洪水 预报 ， 
已 有 的 人 研究 结果 表明 该 模型 具有 较 好 的 模拟 效果 ， 
能 较 准 确 的 估计 发 生 冰 塞 时 的 回 水 剖面 。 本 研究 
选用 了 RIVICE 中 的 水 动力 ` 水 力 糙 率 和 冰 盖 受 力 
增 厚 这 3 个 模块 对 黄河 三 湖 河口 弯 道 处 的 冰 塞 洪水 
过 程 进行 模拟 。 
1.3.1 水 动力 模拟 虽然 河道 上 履 盖 有 冰 ,但 依旧 
将 其 视 为 明渠 水 流 ,忽略 冰 盖 的 压力 。 和 一 般 水 力 
学 模型 相同 ,RIVICE 模型 采用 圣 维 南方 程 组 来 计算 
一 维 明渠 非 恒 定 流 , 因 将 河道 概 化 成 矩形 且 无 侧 向 
入 流 , 圣 维 南 方程 组 被 简化 如 下 : 


连续 方程 BE + 99 -0 (1) 
Ot Ox 


ðt Ox JAg — 


( 
qa A Mpi aL 


运动 方程 : e 十 dr E 
2) 


B ðx 


图 1 黄河 三 湖 河 口 弯 道 段 
Fig. 1 Sanhuhekou bend reach of the Yellow River 
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式 中 :B 为 河道 宽度 (m);z 为 水 位 (m) ;0 为 流量 (mi 
s) ;1 为 时 间 (s) ;x 为 距离 (m) ;ov 为 流速 (ms );4 
为 河道 过 流 面 积 (wm ) ;g 为 重力 加 速度 (9.8 ms 7) ;F 
为 弗 劳 德 数 ;K 为 输送 系数 。4: 为 过 流 面积 随 水 位 
的 变化 率 。 

采用 隐 式 有 限 差分 格式 对 方程 进行 求解 ,得 到 
各 个 计算 断面 上 各 个 时 刻 的 水 流 状 态 。 
1.3.2 水 力 糙 率 在 有 冰 盖 的 河 段 ,水流 除 受 到 河 
床 的 阻力 外 还 受到 冰 盖 下 表面 的 阻力 。 冰 层 履 盖 
下 河道 的 输送 系数 开采 用 调整 后 的 曼 宁 公式 计算 : 

K= LR o) 


WP in WW RAIER A E TEKRAR AVK 
盖 下 的 水 力 半 径 [R=4/(P+P)] ;已 为 河道 湿 周 ; PA 
冰 盖 湿 周 。 复 合 曼 宁 烽 率 系 数 n, 用 Sabaneev 公式 


计算 : 
DERE 2/3 
Zu 3 | (4) 


式 中 :河床 的 曼 宁 糙 率 系数 作为 一 个 参数 ,用 畅 
流 期 的 洪水 进行 率 定 ,在 模拟 冰 塞 的 过 程 中 保持 
不 变 ; 对 于 冰 盖 的 受 宁 烙 率 系数 疡 ,模型 采取 KGS 
TIK”, 

KGS 方 法 基于 Nezhikhovskiy 提 出 的 公式 ,将 
冰 盖 的 米 率 系数 视 为 冰 厚 的 函数 : 
5 esl0.0302 In(h) + 0.0445] (5) 
IAP : h KS ERE BE (m) ; ns a W h=8 m PF YK s Ir] s 
率 系 数 ,将 其 作为 1 个 系数 ,用 冰 塞 时 的 洪水 数据 
1.3.3 冰 盖 厚度 ”作用 在 冰 盖 上 的 力 包 括 6 种 ,如 
图 2 所 示 : 了 为 水 流 对 冰 盖 前 缘 产 生 的 推力 ;为 冰 


m= 


的 分 量 ; Fr 为 冰 盖 和 河岸 之 间 的 摩擦 力 ;Fe 为 当 冰 
盖 冻 结 到 河岸 上 产生 的 内 聚 力 ; 刀 为 冰 盖 的 内 阻 
力 。 因 此 次 模拟 春季 开 河 期 的 冰 塞 ,此 时 严 很 小 ， 


所 以 忽略 。 
水 流 对 冰 盖 前 缘 产 生 的 推力 的 计算 公式 为 : 
Ew Ww 
m= v, Bd, - (6) 


式 中 :4d, 为 冰 盖 前 缘 上 游 处 水 深 (m) ;di 为 冰 盖 前 缘 下 
水 深 (m) ;vw 为 冰 盖 前 缘 下 的 宽度 平均 速度 (ms ');B 
为 冰 羡 宽度 (m);y 为 水 的 比重 (9800 N ^m?) 

冰 羡 下 流动 的 水 对 冰 羡 的 拖 力 的 计算 公 
式 为 : 


dSn}? 
r (05 y (7) 
式 中 :$ 为 水 力 坡度 ;d 为 冰 盖 下 水 深 (m) ;4 为 冰 盖 
的 下 表面 面积 (mw )。 
冰 盖 重量 沿 坡 向 的 分 量 Py 的 计算 公式 为 : 
F,-y,V,S (8) 


式 中 :y 为 冰 的 比重 (9020 N - m?) ;为 两 岸 间 冰 盖 
的 体积 (mi’)。 

冰 盖 和 河岸 之 间 摩 擦 力 fj 的 计算 公式 为 : 

F,-2fhlK, tang (9) 

RP uA Riva B5 FW 73 (Pa) ; Ki 29 17 2828, € 
示 冰 盖 横 向 应 力 与 纵向 应 力 之 比 ;gq 为 冰 盖 和 河岸 
间 的 摩擦 角 ;! 为 断面 之 间 的 距离 (m)。 将 K, 和 tang 
组 合成 1 个 参数 Kitan ,用 冰 塞 时 的 洪水 数据 率 定 。 

冰 羡 内 阻力 也 的 计算 公式 : 


F, i 一 Yer (10) 


式 中 :e 为 冰 盖 的 孔隙 度 ;K; 为 强度 系数 ,二 者 均 作为 


盖 下 流动 的 水 对 冰 盖 的 拖 力 ;4 为 冰 盖 重量 沿 坡 向 


参数 ,用 冰 塞 时 的 洪水 数据 率 定 。 


i a 冰 盖 重力 与 垂直 于 冰 盖 表面 
Bo, 的 重力 分 力 之 间 的 夹 角 


图 2 冰 盖 受 力 示 意图 


Fig. 2 Schematic diagram of ice covert force 
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如 果 +FotFw-Fe<F, 即 内 阻力 可 以 抵抗 外 力 ， 
来 冰 会 在 冰 盖 前 缘 发 生 并 置 之 后 向 上 游 延 伸 ,形成 
REDKE, WR FeFoFFOSFi,BPSEH 7) ASBERIE 
抗 外 力 ,来 冰 发 生 推 挤 使 冰 盖 增 厚 直 到 内 阻力 能 够 
抵抗 外 力 ,形成 冰 坝 。 
14 模型 参数 及 敏感 性 分 析 

参数 的 敏感 性 分 析 是 了 解 参 数 对 模型 影响 大 
小 的 重要 途经 。 对 于 本 次 模拟 来 说 ,研究 冰 塞 水 位 
对 参数 的 敏感 性 ,可 以 揭示 冰 塞 洪水 形成 的 主要 控 
制 因素 ,为 凌 汛 防治 提供 科学 依据 。 

本 研究 选取 传统 的 局 部 敏感 性 分 析 方 法 对 参 
数 和 边界 条 件 进行 敏感 性 分 析 , 基 于 2008 年 .2009 
年 .2010 年 的 数据 ,在 保持 其 他 参数 值 不 变 的 情况 
下 , 单 次 让 1 个 参数 或 边界 条 件 提高 10% ,计算 上 边 
界 处 回 水 水 位 的 变化 。 参 数 敏 感性 用 以 下 公式 
计算 : 


Y- Ypa di 
式 中 :base 表 示 参 数 基 础 值 ;了 表示 状态 变量 , 即 上 
边界 处 的 水 位 ;X 表 示 参 数 ;S, 表 示 变 量 了 对 参数 X 
的 敏感 性 指数 。 


2 结果 与 讨论 


21 畅 流 期 模型 率 定 

首先 ,要 验证 模型 水 动力 部 分 的 有 效 性 以 及 率 
定 参 数 忆 ,研究 选取 2008-08-09 一 08-15 .2009-08- 
22—08-27 这 2 场 洪水 来 率 定 ,再 用 2010-06-10 一 
06-15 的 洪水 数据 进行 验证 。 率 定 和 验证 的 结 
中 , 当 m=0.012 时 ,模拟 结果 较 好 。 图 3 表明 了 2008 
年 .2009 年 和 2010 年 年 3 场 洪水 率 定 和 验证 的 结 
果 , 可 以 看 出 ,模拟 出 的 水 位 与 实测 值 吻合 良好 。 
经 过 统计 计算 ,模拟 的 3 场 洪 水 中 ,水 位 的 确定 性 系 
数 分 别 为 0.91.0.90 和 0.87。 流 量 的 确定 性 系数 均 
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230.99, 

这 个 结果 证 明了 RIVICE 模 型 水 动力 部 分 可 以 
很 好 的 模拟 无 冰 状 态 下 的 洪水 ,为 模拟 冰 塞 洪水 提 
供 了 良好 的 基础 。 
2.2 开 河 期 模型 率 定 

开 河 期 冰 盖 破裂 流 冰 、 卡 冰 以 及 冰 塞 被 冲 向 
下 游 等 过 程 具有 很 大 的 随机 性 ,模拟 如 此 复杂 的 过 
程 造 成 的 回 水 对 于 本 研究 来 说 也 是 巨大 的 挑战 。 
因为 模型 并 不 能 完全 精准 地 捕捉 这 一 复杂 的 变化 
过 程 ,所 以 本 研究 选取 2008 年 .2009 4F. , 2010 EF 
河 期 的 洪水 过 程 来 验证 RIVICE 模 拟 冰 塞 回 水 水 位 
的 可 行 性 。 

开 河 期 流量 的 突然 增 大 会 造成 冰 盖 破裂 和 机 
械 增 厚 ,为 更 好 地 预报 特定 来 冰 条 件 下 一 定 的 流量 
会 造成 的 冰 塞 相应 秦 水 高 度 , 上 边界 条 件 设 为 不 变 
的 入 流 流量 ,下 边界 设置 为 不 变 的 水 位 。 因 为 本 次 
模拟 不 考虑 气温 的 变化 ,所 以 单 次 只 模拟 3 日 的 水 
位 过 程 。 

研究 选择 2008-03-16 一 03-19 .2009-03-16 一 
03-19 .2010-03-22 一 03-25 这 3 场 冰 塞 洪水 事件 对 
RIVICE 模型 进行 率 定 和 验证 ,并 分 别 在 各 年 中 选取 
同 流量 下 畅 流 期 模拟 的 水 位 作为 比较 。 表 1 工 中 列 出 
了 3 次 模拟 的 参数 值 和 边界 值 设 置 。 本 次 模拟 设置 
2 个 边界 条 件 和 9 个 参数 ,边界 条 件 为 上 边界 的 流量 
OO 和 下 边界 的 水 位 丈 。9 个 参数 中 河床 的 曼 宁 糙 率 
系数 凤 `. 计 算 冰 塞 下 表面 曼 宁 糙 率 系数 的 参数 mm 由 
计算 冰 塞 与 河岸 之 间 摩 擦 力 的 参数 Ritan .计算 冰 塞 
内 阻力 的 参数 KK 和 冰 塞 的 孔隙 度 e 在 每 次 模拟 中 保 
寺 不 变 。 冰 塞 前 缘 的 厚度 也 、 冰 塞 下 游 完整 冰 羡 的 


厚度 h 和 每 个 时 刻 冰 的 增 量 Vi 每 年 不 同 ,参考 水 文 
年 鉴 。 此 外 ,模拟 河 段 总 长 为 9.89 km,50 m 设 置 1 
个 计算 断面 , 卡 冰 位 置 x 均 设 为 第 2 个 弯 道 收缩 处 ， 
既 第 167 个 计算 断面 处 。 


2 
时 间 /d 时 间 /d 


图 3 畅 流 期 洪水 水 位 模拟 结果 


Fig.3 Simulation results of flood stage at open water 
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表 1 模型 边界 条 件 和 参数 值 


Tab.1 Boundary conditions and parameter of model 


参数 /边界 条 件 单位 表征 含义 2008 年 2009 年 2010 年 
n, - DES ES 0.012 0.012 0.012 
Msn - ok TERE 0.12 0.12 0.12 
Kitan - 冰 赛 与 河岸 的 摩擦 力 0.218 0.218 0.218 
K, - 冰 塞 的 内 阻力 8.5 8.5 8.5 
e - UKSERS LER EE 0.3 0.3 0.3 
F, m 冰 塞 前 缘 厚 度 0.35 0.25 0.20 
h m 冰 塞 下 游 完整 冰 盖 厚度 0 0 0 
x - 卡 冰 位 置 167 167 167 
V m At" 上 游 来 冰 量 50 60 40 
Q ms 上 游 边 界 流量 800 661 670 
W m 下 游 边界 水 位 1020.7 1020.18 1020.36 


图 4 表明 了 2008 年 .2009 年 .2010 年 开 河 期 模 
拟 的 水 位 过 程 与 实测 的 水 位 过 程 的 比较 以 及 与 同 
流量 下 畅 流 期 模拟 的 水 位 的 比较 。 从 图 中 可 以 看 
出 ,在 开 河 期 ,虽然 模拟 的 水 位 变化 较为 平缓 ,实测 
的 水 位 变化 较为 复杂 ,但 是 模拟 的 水 位 和 实测 值 的 
变化 趋势 一 致 。2009 年 第 2 日 之 后 流量 陡 增 到 
897 ,水 位 过 程 增幅 较 大 ,模拟 的 水 位 不 能 很 好 地 模 
拟 出 这 一 陡 增 的 水 位 。 经 过 统计 计算 , 2008 年、 
2009 年 .2010 年 的 最 大 误差 分 别 为 0.035 m 0.249 
m 0.182 m ,确定 性 系数 为 0.92 .0.79 .0.76。 

模拟 结果 说 明了 模型 可 以 很 好 地 模拟 春季 开 
河 期 由 于 河道 冰 塞 而 造成 的 水 位 升 高 。 与 相同 流 
量 下 畅 流 期 模拟 的 水 位 的 比较 可 知 , 有 冰 塞 影响 时 
的 水 位 比 畅 流 期 的 水 位 高 出 平均 1.4 m, 说 明了 冰 塞 
洪水 的 严重 性 以 及 准确 模拟 并 预报 冰 塞 洪水 的 重 
要 性 。 
2.3 敏感 性 分 析 结 果 

从 敏感 性 分 析 结 果 ( 图 $) 中 可 以 看 出 ,上 边界 
处 的 水 位 对 流量 0 最 为 敏感 ,说 明了 开 河 期 控制 上 


游 来 水 的 重要 性 。 其 次 是 计算 冰 盖 曼 宁 糙 率 系数 
的 参数 ns,、 下 游 边界 水 位 开 。 流 量 的 增 大 除了 会 直 
接 导 致 水 位 升 高 外 还 会 使 并 置 的 冰 盘 发 生机 械 增 
厚 , 冰 羡 的 烟 率 增 大 会 增加 水 流 的 阻力 ,下 游 的 水 
位 升 高 会 提高 歼 水 的 基础 ,这 些 均 是 形成 冰 塞 洪水 
的 重要 成 因 , 所 以 这 些 值 的 增 大 会 造成 冰 塞 回 水 水 
位 的 大 幅 增 大 。 上 游 来 冰 量 的 增加 会 造成 更 加 严 
重 的 冰 塞 ,使 上 游 水 位 雍 高 更 多 。 

计算 冰 塞 内 阻力 的 参数 及 .计算 冰 塞 与 河岸 之 
间 摩 擦 力 的 参数 Kitan 和 冰 塞 的 孔隙 度 e 的 增 大 会 
使 上 边界 处 的 水 位 减 小 ,是 因为 冰 塞 的 内 阻力 和 河 
岸 对 冰 塞 的 摩擦 力 都 是 与 流向 相反 的 力 ,这 2 个 值 
的 升 高 会 增加 冰 塞 的 稳定 性 ,使 冰 塞 的 增 厚 减少 ， 
冰 塞 对 水 流产 生 的 摩擦 阻力 减 小 ,使 回 水 水 位 降 
低 。 在 堵塞 过 程 中 冰 被 压 实 , 孔 际 度 越 大 , 冰 塞 体 
积 减 小 的 越 多 ,造成 的 冰 塞 回 水 水 位 就 越 小 。 

河床 的 烟 率 、 冰 六 前 缘 厚 度 和 冰 塞 下 游 完整 冰 
盖 的 厚度 是 低 敏 感 或 非 敏感 参数 。 

模拟 2008 年 .2009 年 和 2010 年 开 河 期 的 冰 塞 


— 实测 值 ( 开 河 期 ) ---- 模拟 值 ( 开 河 期 ) 一 一 REH) — 河 底 高 程 
1022 1022 1022 
1021 1021 | ®) ns 
g 1020 g 1020 g 1020 
X 1018 36 1018 3Xé 1018 
1017 1017 1017 
1016 mEEGENEE EX. 1016 n 41 —29 5 1016 
时 间 /d 时 间 /d 时 间 /d 


图 4 开 河 期 和 畅 流 期 冰 塞 水 位 模拟 结果 


Fig.4 Simulation results of ice jam stage at break-up period and open water period 
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敏感 度 /% 
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参数 /边界 条 件 
图 5 敏感 性 分 析 结 果 


Fig. 5 Results of sensitivity analysis 


水 位 时 ,依据 实测 数据 设置 的 流量 和 水 位 边界 条 


TE . 冰 盖 厚度 和 上 游 来 冰 量 存在 较 大 差异 ,在 此 影 
响 下 ,参数 和 边界 条 件 的 敏感 性 也 显示 出 较 大 差 
异 , 如 2008 年 和 2010 年 冰 塞 开始 的 位 置 向 下 游 移 
动 ,会 造成 上 游 的 回 水 水 位 升 高 ,而 2009 年 则 相 
反 。 这 充分 说 明了 冰 塞 事件 及 其 造成 的 洪水 水 位 
具有 混沌 性 和 随机 性 。 


jii 


3 结论 


本 研究 用 开源 的 河 冰 模 型 RIVICE 模拟 了 春季 
开 河 期 黄河 三 湖 河口 弯 道 段 的 冰 塞 水 位 ,选取 畅 流 
期 和 开 河 期 的 冰 塞 洪水 对 RIVICE 模型 进行 了 验 
证 ,并 分 析 了 模型 参数 和 边界 条 件 的 敏感 性 ,得 出 
以 下 结论 : 

CD 相同 流量 下 , 相 较 于 无 冰 状 态 , 有 冰 塞 发 生 
会 造成 更 高 的 水 位 ,RIVICE 模型 通过 数值 模拟 的 方 
法 考虑 了 来 冰 在 河道 中 的 卡 堵 . 因 受 力 而 导致 的 机 
械 增 厚 和 冰 对 水 流 的 阻力 作用 ,进而 能 有 效 的 模拟 
出 由 于 冰 塞 导致 的 水 位 抬升 。 

(2) 确定 上 边界 流量 情况 下 ,用 RIVICE 模型 模 
拟 河道 的 冰 塞 回 水 水 位 具有 可 行 性 。 可 以 为 上 游 
水 利 工程 的 调 蓄 提供 科学 依据 ,确保 在 不 会 造成 严 
重 的 下 游 冰 塞 洪水 的 情况 下 合理 下 汇流 量 。 在 没 
水 利 工 程 或 者 水 利 工程 无 法 调节 的 情况 下 ,为 上 游 
来 水 量 造成 冰 塞 洪水 的 严重 性 提供 有 效 的 预警 。 

(3) 冰 塞 造成 的 回 水 水 位 对 流量 最 为 敏感 ,证 
明了 上 游 水 利 工程 对 开 河 期 流量 调控 的 重要 性 ， 
此 水 利 工程 调节 方案 的 有 效 性 是 冰 塞 洪水 防 控 的 
重要 保障 。 
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Numerical simulation of the ice jam stage in the Sanhuhekou 
bend reach of the Yellow River 
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Abstract: In this study, the open-source numerical river ice model RIVICE was used to simulate the ice jam 


water level in the Sanhuhekou bend reach of the Yellow River during the break-up period. Its parameters were 


then subjected to sensitivity analysis. Flow, water level, and ice th ickness were obtained from the Hydrological 
Yearbook. Results show that the RIVICE model can be utilized to simulate the ice jam water level in the break-up 


period. The backwater level of an ice jam is highly sensitive to the upstream boundary discharge, indicating the 


importance of water conservancy projects in flow regulation. Therefore, the established model can provide a 


scientific theoretical basis for managing water conservancy projects and ice flood control. 


Keywords: Sanhuhekou bend reach of the Yellow River; ice jam water level; RIVICE model 


